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and  Aaron  R.  Byrd 


PURPOSE:  This  research  experimentally  determined  phosphorus  mass  distribution  coefficients 
(i.e.,  relationship  between  soil  exchangeable  and  soluble  phosphorus),  exchangeable  phosphorus 
pools,  and  soluble  phosphorus  in  the  interstitial  water  for  soils  exhibiting  a  range  of  phosphorus 
concentrations.  The  research  also  developed  relationships  between  these  variables  and  com¬ 
monly  measured,  crop-available  soil  phosphorus  for  use  in  establishing  initial  soil  P  parameters 
in  the  System-Wide  Water  Resources  Program-Nutrient  Sub-Model  (SWWRP-NSM). 

BACKGROUND:  Mass  transfer  of  soluble  phosphorus  (P)  between  soil  and  laminar  overland 
flow  can  be  described  as: 


sd=K 


c, 


d  2 


(p-c, 


d  y 


(1) 


where 


Sd = the  mass  transfer  flux  (M-L'2-Tx) 
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Cd=  the  concentration  of  soluble  P  in  the  water  column  (M-U  ) 

C/2  =  the  concentration  of  soluble  P  in  the  soil  interstitial  water  (M-U  ;  also  referred  to  as 
the  equilibrium  P  concentration  or  EPC,  Froelich  1988) 
ke  =  the  mass  transfer  coefficient  (L-Tl) 

cp  =  the  porosity  (dimensionless;  Wallach  et  al.  1988;  Chapra  1997;  Gao  et  al.  2004). 

C/2  can  rapidly  (i.e.,  hours)  exchange  with  P  that  is  reversibly  adsorbed  to  soil  constituents 
(primarily  on  Al  and  Fe  hydroxides)  until  an  approximate  equilibrium  is  achieved  between 
particulate  and  aqueous  phases  as: 

C/2<=^2  (2) 

where  Cp2  =  the  soil  adsorbed  inorganic  P  pool  Barrow  1983;  Van  Riemsdijk  et  al. 

1984).  In  general,  Cp2  represents  a  small  fraction  of  the  total  adsorbed  inorganic  P  pool  and  a 
smaller  fraction  of  the  total  potential  capacity  for  soils  to  adsorb  P.  Shifts  in  P  equilibrium  due  to 
crop  uptake  or  fertilizer  application  can  result  in  rapid  (on  the  order  of  minutes  to  hours) 
exchanges  between  the  particulate  and  aqueous  phases.  For  example,  P  desorption  from  soil 
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occurs  when  Cdi  declines  below  equilibrium  conditions  and  P  adsorption  to  soil  occurs  when  Cdi 
exceeds  equilibrium  conditions. 

The  total  concentration  of  this  exchangeable  inorganic  P  pool  {Cn,  M-L~  )  in  the  surface  soil 
layer  can  be  divided  into  two  components  as: 


Ct2  (l  tp)  '  p  '  Cp2  +  (pCd2 


(3) 
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where  p  =  the  soil  density  ( M-L '  ).  Cdi  and  Cp i  are  related  to  an  equilibrium  partition  coefficient 
as: 


Cpl=kd2-Cd2  (4) 

where  kdi  =  the  mass  distribution  coefficient  ( Lr-M  ).  Equations  3  and  4  can  be  combined  as: 

CT2  —  Cd2'[(x  +  <p)  (5) 

where 

a=(\-(p)-p-kd2  (6) 

Alternatively,  Equation  4  can  be  solved  for  Cdi  and  combined  with  Equation  3  as: 

Ct2  —  Cp2  ‘{cc  +  cp)  •  kd 2  (7) 

In  order  to  model  soluble  P  transfer  between  soil  and  overland  flow,  Cn,  Cpi,  Cdi,  and  kdi  need 
to  be  initialized  in  SWWRP-NSM.  This  information  is  not  readily  available  and  usually  must  be 
determined  experimentally  via  Langmuir  isotherm  assays.  Numerous  studies  have  shown  that  P 
mass  transfer  varies  positively  as  a  function  of  the  concentration  of  various  adsorbed  P  pools  and 
the  degree  of  soil  P  saturation  (Pote  et  al.  1996;  McDowell  et  al.  2001;  Torbert  et  al.  2002; 
Vadas  et  al.  2005),  indicating  relationships  between  P  mass  transfer  and  soil  management 
practices.  In  addition,  extraction  techniques  have  been  used  to  approximate  ranges  in  Cdi  and  Cpi 
(see  review  by  McGrechan  2002).  Finally,  Fang  et  al.  (2002)  demonstrated  relationships  between 
the  EPC  and  adsorbed  P  pools.  Little  is  known  about  variations  in  kdi  as  a  function  of  adsorbed 
soil  P  (Lewis  and  McGrechan  2002).  However,  this  information  is  critical  for  soil  mass  transfer 
modeling,  since  Cn,  Cpi,  Cdi,  and  kdi  are  inter-related  (i.e.,  Equations  5  and  7).  In  particular,  kdi 
may  not  be  constant,  but  rather  declines  with  increasing  soil  P  concentration  due  to  decreasing 
availability  of  unoccupied  adsorption  sites.  The  objectives  of  this  research  were  to  1)  estimate 
Cn,  Cpi,  Cdi,  and  kdi  for  soils  with  widely  ranging  P  concentrations,  and  2)  develop  empirical 
relationships  between  these  variables  and  commonly  measured  Bray  and  Mehlich  P  (i.e., 
extractable  crop-available  P)  for  use  in  establishing  initial  soil  P  parameters  for  SWWRP-NSM 
(Johnson  et  al.  2008). 
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METHODS: 

Study  Site.  Eight-Mile  Run  is  a  264-ha  sub-watershed  located  in  the  Upper  Eau  Galle  River 
Basin,  west-central  Wisconsin  (Figure  1).  Livestock  (dairy)  pasture  and  associated  barnyards 
represent  approximately  6  percent  of  the  watershed.  This  land  use  is  located  immediately 
adjacent  to  the  tributary  in  the  lower  portion  of  the  watershed.  Alfalfa  and  com  production 
represent  approximately  49  percent,  and  grass,  CRP,  and  wooded  land  uses  occupy 
approximately  45  percent  of  the  watershed.  Crop  production  occurs  in  the  northern,  eastern,  and 
southeastern  regions  of  the  watershed.  Grass,  CRP,  and  wooded  areas  are  located  in  the  north- 
central  and  western  portions  of  the  watershed.  The  watershed  is  bisected  into  approximately 
equal  areas  by  a  railroad  and  U.S.  Route  12.  Grass,  CRP,  and  wooded  areas  account  for 
61  percent  while  livestock  and  crop  production  represent  69  percent  of  the  land  use  in  the  upper 
and  lower  portions  of  the  watershed,  respectively. 


Figure  1.  A  map  of  the  Eight-Mile  Run  watershed  showing  various  land  uses.  Solid  circles  represent  soil 
sampling  locations. 


Field  Soil  Collection  Procedures.  A  75-m  sampling  grid  was  established  in  the  Eight-Mile 
Run  watershed  for  soil  sampling  between  June  and  early  July,  2006  (Figure  1).  Stations  were 
located  using  high  resolution  digital  orthophotographs  and  differential  GPS  (Garmin  model  72; 
Garmin  International,  Olathe,  KS,  USA).  Vegetation  was  clipped  to  ground  level  at  each 
sampling  point  and  six  5-cm-deep  soil  cores  were  collected  around  a  2-m-diam  circle.  The  soil 
cores  collected  at  each  station  were  combined  into  one  composited  sample.  Land  use  was 
recorded  at  the  time  of  sampling.  In  the  laboratory,  soils  were  passed  through  a  2-mtn  mesh 
screen  to  remove  root  material  and  air  dried  at  25  °C.  Overall,  approximately  430  composited 
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soil  samples  were  collected  from  the  Eight-Mile  Run  watershed.  The  dried  soil  samples  were 
further  composited  for  chemical  analyses  as  a  function  of  37  fields  exhibiting  homogeneous  land 
uses  (Figure  1). 

Soil  Phosphorus  Pools  and  Equilibrium  Characteristics.  Soil  was  extracted  in 
0.025  M  HC1  and  0.03  M  NH4F  and  a  solution  containing  0.2  M  CH3COOH,  0.25  M  NH4NO3, 
0.015  M  NH4F,  0.013  M  HNO3,  and  0.001  M  EDTA  for  determination  of  crop-available  Bray 
and  Mehlich-3  P,  respectively  (Pierzynski  2000).  Langmuir  isotherm  assays  were  conducted 
according  to  methods  described  in  Nair  et  al.  (1984).  Soils  were  subjected  to  initial  P  (as 
KEEPO/t-P  in  a  0.01  M  CaCE  solution)  concentrations  of  0,  0.125,  0.25,  0.50,  1,  5,  10,  25,  50  and 
100  mg-L'1  using  a  soiksolution  ratio  of  1:25  (40  g-L'1).  The  soil  solution  tubes  were  gently 
shaken  in  a  darkened  environment  at  20  °C  over  a  24-hr  period.  The  equilibrated  samples  were 
centrifuged  at  500  g  and  filtered  through  a  0.45-p  filter  for  soluble  reactive  P  (SRP) 
determination  (American  Public  Health  Association  (APHA)  1998). 

The  change  in  SRP  mass  [(initial  Cd2  -  final  Cd2)'Solution  volume]  was  divided  by  soil  mass  to 
determine  the  mass  of  P  desorbed  or  adsorbed  from  soil  (ACP2;  mg-kg'1).  ACP2  was  plotted  as  a 
function  of  final  Cd2  (he.,  equivalent  to  Cd2  in  Equation  1)  to  determine  the  mass  distribution 
coefficient  (kd2;  L-kg'1),  the  equilibrium  phosphorus  concentration  (EPC,  the  point  where  net 
sorption  is  zero  and  equivalent  to  Cd2;  Froelich  1988),  and  the  concentration  of  the  initial 
adsorbed  P  pool  (Cp;  mg-kg'1).  kd2,  Cd2,  and  CP2  (i.e.,  equivalent  to  CP2  in  Equation  1)  were 
calculated  via  regression  analysis  (Statistical  Analysis  System  1994)  from  linear  relationships 
between  final  Cd2  and  ACP2  near  the  point  where  net  sorption  was  zero  (i.e.,  crossover  point; 
Figure  2).  kd2  was  equal  to  the  slope  of  the  regression  line  and  CP2  was  estimated  as  the  intercept. 
This  calculation  was  modified  for  soils  with  very  high  P  content  because  the  ACp2  versus  Cd2 
relationship  was  nonlinear  near  the  crossover  region  and  linear  regression  of  these  data  would 
have  resulted  in  an  underestimate  of  CP2  (Figure  3).  Under  these  conditions,  Cp2  was  estimated 
via  regression  analysis  at  very  low  Cd2  and  kd2  was  calculated  as  the  slope  of  the  line  between  CP2 
and  positive  ACP2  (i.e.,  net  adsorption  of  P)  nearest  the  crossover  point. 

The  P  sorption  capacity  (PSC;  M-L'  )  was  estimated  as  ACp2  for  soils  subjected  to  a  100  mg-L' 
solution  (Nair  et  al.  1998).  The  degree  of  P  saturation  index  (DPSindex;  %)  was  calculated  as: 


DPSIndex  =  100. 


Adsorbed 


PSC  +  P 


Adsorbed 


(8) 


where  P Adsorbed  was  either  Bray  or  Mehlich  P. 
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b) 

6) 

E 


cd2  (mg-L-1) 


Figure  2.  An  example  of  an  adsorption-desorption  isotherm  ■  Cd2  represents  the  concentration  of  soluble 
reactive  phosphorus  in  the  aqueous  phase  after  continuous  shaking  with  soil  over  a  24-hr 
period.  ACp2  represents  the  change  in  soluble  reactive  phosphorus  concentration  after  the 
24-hr  shaking  period  (initial  Cd2  -  final  Cd2)  normalized  with  respect  to  soil  mass.  A  negative 
ACp2  indicates  desorption  of  phosphorus  from  soil  to  the  aqueous  phase  while  a  positive  ACp2 
represents  adsorption  of  aqueous  phosphorus  onto  soil  particles.  The  crossover  point  is 
equivalent  to  the  equilibrium  phosphorus  concentration  (EPC).  The  slope  of  the  regression 
equation  near  the  crossover  point  is  the  mass  distribution  coefficient  (kd2;  Lkg"1).  The  intercept 
of  the  regression  equation  represents  initial  adsorbed  soil  phosphorus  concentration  (Cp2; 
mg-kg"1). 
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Figure  3.  Adjustment  used  to  determine  kd2,  Cd2,  and  Cp2  for  soils  with  very  high  soil  Bray  and  Mehlich-3 
P  concentrations  and  nonlinear  patterns  near  the  crossover  point.  Cp2  is  estimated  as  the 
intercept  for  the  linear  portion  of  the  isotherm  at  low  Cd2  concentrations.  kd2  is  calculated  as 
the  slope  of  a  line  passing  through  Cp2  and  the  first  data  point  above  the  crossover  (i.e., 
positive  ACp2).  The  red  line  represents  the  adjusted  kd2  slope. 


RESULTS  AND  DISCUSSION:  Bray  and  Mehlich  P  exceeded  30  mg-kg'1  and  50  mg-kg'1 
(optimum  for  crop  growth),  respectively,  primarily  in  the  lower  portion  of  the  watershed  in 
conjunction  with  crop  and  livestock  production  land  use  practices  (Figure  4).  CRP,  grasses,  and 
woodlot  land  uses  were  associated  with  Bray  and  Mehlich  P  concentrations  that  were  at  or  below 
optimum  levels  for  crop  uptake. 

Examples  of  ACP2  versus  Cd2  are  shown  for  four  fields  in  Figure  5.  High  kd2  was  associated  with 
low  Cd2,  crop-available  P,  and  CP2.  It  declined  in  conjunction  with  higher  soil  P  concentrations. 
For  all  fields,  kd2  covaried  nonlinearly  as  a  function  of  crop-available  P  (Figures  6a  and  6b).  It 
was  greatest  for  soils  with  low  Bray  and  Mehlich  P  and  declined  in  a  negative  logarithmic 
pattern  with  increasing  soil  P.  These  patterns  may  be  related  to  differences  in  the  degree  of  P 
saturation  of  sorption  sites  (Kleinman  et  al.  2000).  For  instance,  the  soil  buffering  capacity  to 
adsorb  P  under  conditions  of  P  disequilibrium  would  decline  as  DPS  increases  due  to  a  lower 
number  of  unoccupied  adsorption  sites,  resulting  in  lower  P  sorption  efficiency  (see  below  and 
Figure  7). 
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Figure  4.  Spatial  variations  in  soil  Bray  and  Mehlich-3  P  concentrations  in  Eight-Mile  Run  in  2006. 
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kd2=  40  L-kg-1; 


Cd2  =  0.013  mg-L1;  Cp2  =  7.3  mg-kg1;  Bray  P  =  15  mg-kg1;  Mehlich  P  =  22  mg-kg1 

Cd2  =  0.367  mg-L1;  Cp2  =  15.4  mg-kg1;  Bray  P  =  79  mg-kg1;  Mehlich  P  =  78  mg-kg1 

Cd2  =  0.990  mg-L1;  Cp2  =  53.7  mg-kg1;  Bray  P  =  290  mg-kg1;  Mehlich  P  =  286  mg-kg1 

Cd2  =  2.688  mg-L-1  ;CP2  =  106.3  mg-kg1;  Bray  P  =  350  mg-kg1;  Mehlich  P  =  382  mg-kg1 
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Figure  5.  Variations  in  the  soil  mass  distribution  coefficient  (kd2),  equilibrium  soluble 
phosphorus  concentration  (Cd2),  and  soil  adsorbed  phosphorus 
concentration  (Cp2)  for  four  soils  exhibiting  different  Bray  and  Mehlich-3 
phosphorus  concentrations. 
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Figure  6.  The  soil  mass  distribution  coefficient  (kd2;  panels  a  and  b)  and  the  equilibrium  soluble 

phosphorus  concentration  (Cd2;  panels  c  and  d)  versus  soil  Bray  and  Mehlich-3  phosphorus. 
The  NLIN  procedure  (Statistical  Analysis  System  1994)  was  used  to  estimate  the  coefficients 
a,  bi,  b2,  and  X0.  See  Equations  7  and  8  for  a  description. 


Figure  7.  Initial  adsorbed  soil  phosphorus  (Cp2)  versus  soil  Bray  and  Mehlich-3  phosphorus. 
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Cd2  exhibited  a  biphasic  linear  increase  as  a  function  of  crop-available  P  (Figures  6c  and  6d). 
These  patterns  could  be  interpreted  using  a  segmented  linear-linear  regression  analysis  (NLIN; 
Statistical  Analysis  System  1994)  that  estimates  a  threshold  between  regression  lines  (Kleinman 
et  al.  2000).  Linear  models  are: 

y  =  a  +  {bx-x)  (9) 


y  =  a  +  (brx0)  +  b2-(x-x0 ) 


(10) 


where  a,  bi,  and  b2  are  constants  and  x0  is  the  Bray  or  Mehlich  P  threshold  concentration.  x0  was 
similar  for  Bray  and  Mehlich  P  at  183.7  and  180.2  mg-kg'1,  respectively.  At  these  threshold 
levels,  Cd2  was  ~0.4  mg-L'1  for  Bray  and  Mehlich  P,  respectively.  Above  these  levels,  the  rate  of 
change  in  Cd2  with  respect  to  Bray  or  Mehlich  P  increased  (Figures  6c  and  6d).  Similar  to  kd2, 
these  patterns  are  likely  related  to  increasing  DPS.  As  sorption  sites  become  increasingly 
saturated,  buffering  capacity  declines,  resulting  in  greater  equilibrium  Cd2-  Strong  linear 
relationships  were  observed  between  Bray  or  Mehlich  P  and  CP2  (Figure  7).  Overall,  CP2 
represented  approximately  25  percent  of  the  Bray  and  Mehlich  P  concentration. 


Land  uses  associated  with  crop  and  livestock  production  in  the  lower  portion  of  the  watershed 
exhibited  the  highest  DPSBray  and  DPSMehiich  (Figure  8).  Both  indices  approached  80-percent 
saturation  in  conjunction  with  high  Bray  and  Mehlich  P  (Figure  4).  kd2,  Cd2,  and  CP2  covaried 
similarly  with  increasing  DPSBray  and  DPSMehiich  (Figure  9).  kd2  declined  logarithmically  while 
Cd2  increased  in  a  linear  biphasic  pattern  with  increasing  DPSBray  and  DPSMehiich-  Threshold 
values  for  DPSBray  and  DPSMehiich  were  37.3  percent  and  39.1  percent,  respectively,  coinciding 
with  a  Cd2  of  -  0.4  mg-L'1.  Regression  relationships  between  CP2  and  DPSBray  or  DPSMehiich  were 
linear  and  positive  (Figure  9). 
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Figure  8.  Spatial  variations  in  the  degree  of  phosphorus  saturation  (DSP),  based  on  soil  Bray  and 
Mehlich-3  phosphorus  concentration. 
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Figure  9.  The  soil  mass  distribution  coefficient  (kd2;  panels  a  and  b),  equilibrium  phosphorus 

concentration  (Cd2;  panels  c  and  d),  and  the  initial  adsorbed  soil  phosphorus  concentration 
(Cp2;  panels  e  and  f)  versus  the  degree  of  phosphorus  saturation  (DPS)  based  on  soil  Bray 
and  Mehlich-3  phosphorus  concentration. 
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An  important  finding  of  this  study  is  that  commonly  measured  crop-available  P  can  be  used  to 
initialize  CP2,  Cd2,  and  kd2  in  SWWRP-NSM.  Empirical  regression  equations  are  summarized  in 
Table  1.  Relationships  between  observed  and  predicted  variables  were  highly  significant  and 
slopes  were  near  1.0  (Figure  10).  Ct2  can  be  estimated  from  Equation  3.  Additional  information 
on  Cp2,  Cd2,  and  kd2  is  needed  for  a  variety  of  soil  types  in  order  to  develop  general  empirical 
patterns  for  initializing  these  variables  in  SWWRP-NSM. 


Table  1 

Empirical  Equations  for  Initializing  kd2,  Cd2,  and  Cp2  in  SWWRP-NSM  Using 
Commonly  Measured  Soil  Bray  or  Mehlich-3  Phosphorus  (P). 
x  =  Bray  or  Mehlich  P  (mg/g) 

Dependent  variable,  y 

Bray  P  regression  equations 

Mehlich  P  regression  equations 

kd 2,  L  kg  1 

ln(y)  = -0.661 -ln(x)  + 7.960 

For  x  <  183.7  mg-kg'1  Bray  P 
-  0.02067  +  0.00208-x 

ln(y)  =  -0.850-ln(x)  +  8.750 

For  x  <  180.2  mg-kg  1  Mehlich  P 
-  0.04982  +  0.00267-x 

Cd 2.  mg-L'1 

Forx  >  183.7  mg-kg  1  Bray  P 

0.36118  +  0.0080-(x-  183.7) 

For  x  >  1 80.2  mg-kg'1  Mehlich  P 

0.4297  +  0.0064-(x-  180.2) 

Cp2,  mg-kg  1 

y  =  0.252-x  +  1.627 

y  =  0.223-x  +  3.939 

SUMMARY:  Langmuir-type  isothenn  soil  assays  were  conducted  to  estimate  the  mass 
distribution  coefficient  (kd2),  the  soil  adsorbed  inorganic  P  pool  (CP2),  and  the  soluble  P  pool  in 
the  soil  interstitial  water  (Cd2)  over  a  range  of  soil  P  concentrations  for  use  in  initializing  soil  P 
compartments  in  SWWRP-NSM.  kd2  was  not  constant;  but  rather,  varied  in  a  negative 
logarithmic  pattern  as  a  function  of  increasing  crop-available  P  concentration.  This  pattern  may 
be  attributable  to  declining  soil  buffering  capacity  for  P  as  DPS  increases  and  availability  of 
unoccupied  adsorption  sites  decreases.  CP2  increased  in  a  linear  pattern,  while  Cd2  exhibited  a 
biphasic  linear  increase,  in  conjunction  with  increasing  crop-available  P.  The  latter  pattern  may 
also  be  related  to  the  DPS. 
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